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D ■ Abstract. This article is devoted to computing the lower and upper bounds of the Laplace 



eigenvalue problem. By using the special nonconforming finite elements, i.e., enriched 
Crouzeix-Raviart element and extension Q\ ot , we get the lower bound of the eigenvalue. 
Additionally, we also use conforming finite elements to do the postprocessing to get the 
upper bound of the eigenvalue. The postprocessing method need only to solve the corre- 
sponding source problems and a small eigenvalue problem if higher order postprocessing 
method is implemented. Thus, we can obtain the lower and upper bounds of the eigenvalues 
simultaneously by solving eigenvalue problem only once. Some numerical results are also 
presented to validate our theoretical analysis. 



■ 1. Introduction 

t> 

Q\ ■ The eigenvalue problems are important, which appears in many fields, such as quantum 
^ ■ mechanics, fluid mechanics, stochastic process and etc. Thus, a fundamental work is to find 
the eigenvalues of partial differential equations. From last century, abundant works are 
dedicated to this topic. 

Feng in his famous paper [6] cites the pioneer work of Polya in computing the upper bound 
of Laplace eigenvalue problem. And based on the minimum-maximal principle discovered 
by Rayleigh, Poincare, Courant and Fischer ect., any conforming finite element method will 
give the upper bound (see Strang and Fix [IB])- Nevertheless, to the lower bound aspect, 
until 1979, Rannacher [T7] gives some numerical results for plate problem. 

And then there is few work on analysis of the lower bound for a long time. Hu, Huang and 
Shen [7] get the lower bound of Laplace equation by conforming linear and bilinear elements 
together with the mass lumping method. Inspired by the minimum-maximal principle, 
people try to find the lower bound with the nonconforming element methods. Recently, 
a series of works make progress in this aspect, e.g. Lin and Lin [H] use the asymptotic 
expansion skill to compute the eigenvalues by nonconforming finite element method; also 
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see the numerical reports of Liu and Liu pj)], Liu and Yan [16], and the work of Lin, 
Huang and Li [10]. Another way by Armentano and Duran [Tj is to use another kind of 
expansion method to get the lower bound, which is of less restriction to the smoothness of 
eigenfunctions compared with the asymptotic expansion skill. Also see the follow-up works 
by Li [9], Lin p], Yang [21]. 

Inspired by these works, this article propose a method to obtain the lower and upper 
bounds of the eigenvalue simultaneously which only need to solve the eigenvalue problem 
once and additional auxiliary source problem. Our method can be described as follows: (1) 
solve the eigenvalue problem by some nonconforming finite element; (2) solve an additional 
auxiliary source problem in an conforming finite element space. Since we obtain not only 
lower bound but also upper bound of the eigenvalue, we can give the accurate error of the 
eigenvalue approximations. Compared with the existed literature, this note contributes on 
the following aspects: 

• The assumption of the lower bound in [1] needs an critical assumption that: \u — 
Uh\i > Ch 1 , 7 < 1, which promise itself be the dominant term in the expansion. By 
our recent result of lower bound of convergence rate by finite element method, we 
get rid of this constrain. 

• A new application of the correction method of eigenvalue problem is proposed to 
obtain the lower and upper bounds of the eigenvalues by solving the eigenvalue 
problem once and an additional source problem. 

• By using higher order conforming finite element to do the correction in a new way, 
we also prove the upper bound of the corrected eigenvalue approximations. 

• After calculating upper and lower bounds simultaneously, we can find how much 
accuracy we have actually achieved (an accurate a posteriori error estimate). 

For simplicity, we only discuss the problem in TZ 2 , but the methods and results here can be 
easily extended to the case TZ 3 . In this paper, we will use the standard notation of Sobolev 
spaces ([IE|). The outline of the paper is as follows. In Section 2, some preliminaries 
and notation are introduced. The weak form of the Laplace eigenvalue problem and its 
corresponding discrete form is stated. In Section 3, we will give the results about the 
lower bound with nonconforming finite elements. Section 4 is devoted to analyzing the 
upper bound of the eigenvalue by postprocessing with the lowest order conforming finite 
element methods. In section 5, another type of postprocessing method is proposed to obtain 
not only higher order accuracy but also upper bound approximation of the eigenvalues. 
Some numerical results are presented in section 6 to test our theoretical results and some 
concluding remarks are given in the last section. 

2. The eigenvalue problem 
In this paper we are concerned with the Laplace eigenvalue problem: 
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Find (A, u) such that 

{—An = An, in tt, 
u = 0, onffi, 
J n u 2 dxdy = 1, 

where Q is a bounded domain in 1Z 2 with continuous Lipschitz boundary dQ . 
The variational problem associated with (12. II) is given by: 
Find u G V := H^(Q) and A G TZ such that b(u, u) = 1 and 

(2.2) a(u,n) = A6(«,n), Vi> G V, 
where 

a(n,t>) = / Vu-Vvdxdy and b(u,v) = / uvdxdy. 
Jn Jn 

From [3], we know that the Lapalce eigenvalue problem has a positive eigenvalue sequence 
{Aj} with 

< Ai < A 2 < ■ ■ • < Afc < • ■ • , lim X k = oo, 
and the corresponding eigenfunction sequence {uj} 

Ui,u 2 , ■ ■ ■ ,u k , ■ ■ ■ , 

with the property b(ui,Uj) = 5ij. 

From the result in [3], the Rayleigh quotient is defined by 

t-> / \ a(u,u) 

RW := = A. 

b{u, u) 

We define the discrete finite element method to solve the problem (12. 3p . Th is a quasi- 
uniform triangulation. Based on this partition, Sh denotes the set of all edges in partition 
Th- The finite element space Vh is the corresponding finite element space to the partition, 
i.e. V h N ^Fasa nonconforming space and V h C V a conforming space. 

In the rest of this paper, the finite element space Vh can be V h c or V h NC . The finite element 
approximation of (12.21) is defined as follows: 

Find (Xh,Uh) ElZxVh such that b(uh,Uh) = 1 and 

(2.3) a h (u h ,v h ) = X h b(u h ,v h ), Wv h e V h , 

where the bilinear form a^(-, •) coincides with a(-, ■) in conforming finite element, or a ele- 
mentwise representation of a(-, ■) in nonconforming situation, e.g. 



ah{Uh,v h ) = 

KdT h 



^ / Vu h ■ Wv h dxdy. 

J k 
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Obviously, the bilinear form for conforming situation can also be presented as this form. It 
is easy to see that for both situation, the bilinear are Vh— elliptic. Thus, we define the norm 
on 14 + V by 



IL = a h (v,v), forveV^ c , 



and 



\v\\l = a(v,v), ioiveV + Vh- 



For both conforming and nonconforming situations, the Rayleigh quotient holds for the 
eigenvalue A^ 

a h (u h ,u h ) 

R \ u h) = T7 T = X h- 

b{u h ,u h ) 

Similarly, the discrete eigenvalues problem ( 12. 3 p has also an eigenvalue sequence {Xj,h} with 
and the corresponding discrete eigenfunction sequence {uk,h} 



u l,h, U2,hi " " " j u k,h, ■ ■ ■ , UN,h 

with the property b(v,i^, = 1 < i, j < N h (N h is the dimension of 14)- 

We cite here the classical result about eigenvalue: minimum-maximum principle. Let 
Aj be the j— th eigenvalue of (12. 2ft and A^ be the j— th eigenvalue of (12. 3p . respectively. 
Arranging them by increasing order, then we have ([3]) 

(2.4) A,- = min maxR(v), X~h = mm maxR(v). 

V,cy,dimV,=i veVj J ' VjCV h ,<iimVj=j veVj 

Since the convergence of the finite element approximation to the eigenvalue problem 
depends on the regularity of the original eigenvalue problem, here we assume that the 
regularity of the eigenfunction u G H 1+ "'(Q) with < 7 < 1 which is decided by the largest 
inner angle of the boundary dfl. 

For the eigenvalue problem, we have the following basic expansion from pQ , and has been 
extensively used in [131 EH [2T] . 

Lemma 2.1. Suppose (X,u) is the eigenpair of the original problem (OOP, (Xh, u^) G 71x14 
is the eigenpair of the discrete problem A2.3\) , we have the following expansion 

X-X h =\\u- u h \\l fh - X h \\v h - u h \\l 

(2.5) +X h (\\v h \\ 2 b - \\u\\ 2 b ) + 2a h (u- v h ,u h ), Mv h eV h . 
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3. Lower bound with nonconform finite element methods 

In this paper, we are concerned with two types of nonconforming finite elements: Enriched 
Crouzeix-Raviart (ECR) ([3 [13]) and Extension Q[ ot (EQ{ ot ) ([12]) for triangle and rectangle 
partitions, respectively. 

• ECR element is defined on the triangle partition and 

V h C '-={ v e L 2 {Q) : v\ K g span{l,x,y,x 2 + y 2 }, J v\ Kl ds = Jv\ K2 ds, 

(3.1) when K 1 nK 2 = £e £ h , and J v\ K ds = 0, if K n dQ = ij, 

where K, Ki, K 2 G %■ 

• On the rectangle partition, EQ 1 ^ is employed defined by 

v h° '-={y e L 2 (Q) :v\ K e span{l, x, y, x 2 , y 2 }, J v\ Kl ds = J^v\ K2 ds, 

(3.2) if K x n K 2 = £, and J v\ K ds = 0, if K n dtt = ij, 

where K, K\, K 2 € %■ 

From [8] and [H], the following basic error estimates for the two nonconforming finite 
elements hold 

(3.3) \X-X h \ < Ch^\\uf 1+1 , 

(3.4) \\u-u h \\ aih < C/i 7 ||u|| 1+7 , 

(3.5) \\u-u h \\ h < Ch?\\u- u h \\ a)h <Ch 2l \\u\\ 1+1 . 

The interpolation operator corresponding to ECR and EQ^ can be defined in the same 
way: 

(3.6) i(u-n h u)ds = 0, Weffc, 

J i 

(3.7) f (u- n h u)dK = 0, WK G 7;, 

Jk 

Lemma 3.1. (See [X3\) For the ECR and EQ™* elements, the interpolation operator satis- 
fies: 

\\u - II h u\\ b + h\\u - n h u\\ a ^ h < C7i 1+7 ||w||i +7 , 

for any u G H 1+1 (Q). 

Lemma 3.2. With the interpolation defined above, for any u G V , we have the following 
results 

(3.8) a h (u - II h u, n h u) = 0, 

(3.9) \\n h u\\ ath < C\\u\\ a , 
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(3.10) ||u - II h u\\ b < Ch\\u - n h u\\ a ,h, 

(3.11) \\u - II h u\\ b < Ch\\u - U h \\ a ,h- 

Proof. Here we only give the proof for the ECR element and the one for EQ^ is almost 
the same. 

Integrating by parts, we have 

a h (u - II h u, n h u) = [ / (u - n h u)AII h udxdy + / (u - n h u)(Vn h u) ■ nds]. 
KeT h Jk JdK 

Since U h u G V^ G 

AllhU = const. 

From the definition of the face interpolation, the following equality holds 

(u — n h u)AII h udxdy = 0. 



K 



For any £ G OK can denoted by the linear function y = kx + m, then the normal direction 
Dr 

y 2 ), we have 



corresponding to this edge is n = -j==(±k, —1). Thus suppose IJ h u = a + bx + cy + d(x 2 



(Vn h u)-n = ±(b + 2dx,c + 2dy)^^=(k,-l) rj 

yl + k 2 

(bk — c — 2dm) = const. 



VT+W 

From the definition of the interpolation, we have 



(u - n h u){Vn h u) ■ nds = 0. 

J i 



So 

a h (u - IJ h u, n h u) = 0. 

From (I3.8p . we can obtain 

a h (n h u,n h u) = a h (u,n h u) < \\n h u\\ a)h \\u\\ ath . 

Canceling the term ||i7^u|| a ^ leads to ( 13.91) . 

Defining il be the piecewise constant interpolation, we have 

\\u-n h u\\ 2 b = b(u - n h u,u - n h u - n (u - n h u)) 

< \\u - n h u\\ b \\u - n h u - n (u - n h u)\\ b 

< c\\u - n h u\\ b h\\u - n h u\\ a ,h- 

It means we arrive the result (13.101) . 

For the last inequality, from (13.91) and (13.101) . we have following equality 

\\u - IJ h u\\ b < Ch\\u - n h u\\ a ,h < Ch(\\u - u h \\ a , h + \\u h - II h u\\ a , h ) 
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= C7i(||u - u fc || 0)ft + \\U h u h - n h u\\ a>h ) < Ch\\u - u h \\ a>h . 
Thus, we get fl3TTD . □ 

Lemma 3.3. Section 3}) If we solve the eigenvalue problem Ii2.2\) by , ECR, bilinear 
or linear elements, the following lower bound for the convergence rate holds 

(3.12) \\u-Uh\\a,h > Ch. 

Theorem 3.4. Let Xj and be the j-th exact eigenvalue and its corresponding numerical 
approximation by ECR or EQ r ° l element. Assume uj e H 1+1 (Q) with < 7 < 1. When h 
is small enough, we have 

(3.13) < Xj - X jth < Ch 2j \\uf 1+r 

Proof. The result for ECR and EQ^ elements can be proved in the uniform way. We choose 
Vj t h = IlhUj in Lemma [2.11 For the second term in (12. 5ft . from Lemma [3.21 we have 

\\uj,h - n h Uj\\l < Ch 2 \\uj -Uj jh \\l th - 
For the third term in (I2.5|) . from (13.111) . we have 

\\ V 3,h\\b ~ \\ U j\\b = (n h Uj -Uj,n h Uj + Uj) 

= {nhUj - uj, (n h u jth + uj) - n (n h u jih + Uj)) 

< Ch\\n h Uj - Uj\\ b < Ch 2 \\uj - Uj th \\ a>h , 

where Ilo denotes the piecewise constant interpolation. Together with (13 .8p and Lemma l3~3l 
the first term will be the dominant term in ( 12.5ft and then ( 13.131) can be derived. □ 

4. Upper bound with linear conforming element 

4.1. Upper bound by direct use of linear conforming element. Due to the minimum- 
maximal principle (12.41) . by using conforming finite element, we can get the upper bound 
naturally. Here we use the lowest order conforming element for both triangle and rectangle 
partitions, respectively 

• On the triangle partition, we use the linear element defined as below: 

(4.1) Vf := {t; G C°{Q) : v\ K G span{l,x,y}} DH^tl). 

• On the rectangle partition, the bilinear element is employed here: 

(4.2) Vf:= [veC°(f2) :v\ K es P an{l,x,y,xy}}nHi(£l). 

In the conforming finite element space , we can define the corresponding discrete 
eigenvalue problem: 

Find (Xh, Uh)€ TZ x Vff such that b(uh,Uh) = 1 an d 

(4.3) a(u h ,v h ) = X h (u h ,v h ), \/v h E '. 



(4.7) 


< A - X h 


< 


(4.8) 


< X h - A 


< 


(4.9) 


o < x h - x h 


< 
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From [3], the basic estimates about approximation errors exist 

(4.4) 0<A,-A < Ch^\\u\\\ + ^ 

(4.5) \\u~ U h \\ a , h < C7l 7 ||tt||i +7 , 

(4.6) \\u-Uh\\ b < Ch^\\u\\ 1+7 . 

So by (14. 4p . (I3.3P and Theorem 13. 4[ we have the following estimate between conforming 
and nonconforming eigenvalue approximations: 

Ch^\\u\\\ +V 
Ch^\\u\\l + , n 
Ch^\\u\\l +r 

4.2. Upper bound by postprocess. The results (14. 71) - (14. 9 p are very interesting and use- 
ful. Especially, (14. 9 p can give a guaranteed error for the current numerical approximations. 
Unfortunately, in order to get (I4.9p . we need to solve the eigenvalue problem twice, one 
with nonconforming and the other with conforming elements, which is always difficult since 
solving eigenvalue problem need much more computation than solving the corresponding 
source problems. So, in this subsection we give a postprocessing method to obtain the upper 
bound eigenvalue approximation without solving eigenvalue problem but only need to solve 
a source problem. 

After obtaining the eigenpair approximation (Xh,Uh) E 7Z x V h by the nonconforming 
finite element, we put them as the right hand side of an auxiliary source Laplace problem. 
Then we utilize conforming finite element method to solve this problem which is defined as 
follows: 

(4.10) a(u h ,v h ) = X h b(u h ,v h ), Wv h E V h c 

After obtaining u h) we calculate the following Rayleigh quotient as an approximation of A 

(4.11) X h :=R(u h ). 

In order to analyze the error estimate for the eigenpair approximation (Xh,Uh) E 71 x V h c , 
we define the projection operator associated with the space V h c as follows 

(4.12) a(P h u,v h ) = a(u,v h ), Wv h E V h c . 
For this projection operator, we have the following error estimate 



(4.13) \\P h u-u\\ b + W\\P h u-u\\ a < Ch 2 ~<\ 



it 



1+7- 



An important lemma from Babuska and Osborn ([3j [19]), tells that 
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Lemma 4.1. (^[3j Lemma 4.1],) For the self-adjoint problem ^2. 3\) . suppose (X,u) be the 
exact eigenpair. Then for any w G V , \\w\\b ^ 0, the Rayleigh quotient R(w) satisfy 

1 2 II 1 1 2 

(4.14) R(w) - A = "~ 7 a - A 



w — u\\„ , \\w — U\\ b 



>\\l \H\b 



Theorem 4.2. For the eigenpair approximation (Aft,«ft) G 1Z x V^p, i/ie following error 
estimates hold 

(4.15) ||u-«ft||o < C^ 7 ||^||i+7, 

(4.16) | A h — A| < Ch 2 \ 

Assume u G H 1+1 {VL) with < 7 < 1. Then when h is small enough, we have the following 
upper bound for Xh defined by U-19{ ), 

(4.17) X h > X. 
Proof. First, we have the following estimate 

\\u h - P h u\\l = a(u h - P h u, u h - P h u) = b(X h u h - An, u h - P h u) 

< \\X h U h - Xu\\ Qjh \\u h - P h u\\ a 

< (\X h \\\u - U h \\ b + I Aft - A||H| b ) \\u h - P h u\\ a 

< Ch 2l \\u h -P h u\\ a \\u\\ 1+1 . 
So the following estimate holds 



(4.18) \\u h -P h u\\ a < Ch 2 ^\\ 



u 



1+7- 



Then we have the following error estimates for 

\\u h -u\\ a < \\u h - P h u\\ a + \\P h U - U\\ a 

< CW\\u\\ 1+ „ 

\\u h -u\\b < \\u h - PhU\\ a + \\PhU ~ U\\ b 

< Ch 2 ^\\u\\ 1+T 
Since V h C V, replacing w with Uh in ( 14. 14j) . we have 

(4-19) A,-A = ^^-A^^. 

Applying the estimates of Uh to (14.191) . we can obtain (14.161) . Furthermore from Lemma [331 
we have \\u — Uh\\ a > Ch. Thus we can see that the first term is dominate and the desired 
result (I4.17P is derived. □ 

Theorem 4.3. Under the conditions in Theorem 3.J\ and \4-2\ for the eigenvalue approxima- 
tion Aft and Xh, we have < X h — X h < Ch 21 and max{|A — X h \AX h — X\} < X h — X h < Ch 2j . 
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Proof. Based on the results in Theorems 13.41 and 14.21 and the error estimate f l 3 . 3 H . we can 
easily prove this theorem. □ 



5. Better upper bound approximation with finer finite element space 

In the last section, the postprocessing method is applied to get the upper bound of the 
eigenvalue which has the same convergence order as the lower bound. But as we know from 
[I~9] . we can employ the postprocessing method to obtain a new eigenpair approximation 
with better accuracy than the obtained nonconforming approximation (A^, Uh). But because 
of the higher order convergence, the analysis of upper bound in last section can not be used 
in this case. Here, we propose a method to produce not only higher order accuracy but 
also upper bound approximation of the eigenvalue. This procedure contains solving some 
auxiliary source problems and a very small eigenvalue problem. 

The aim of this section is to obtain the approximation of the first m eigenvalues A x < 
A 2 < • ■ ■ < A m . Assume we have obtained the first eigenvalue approximations X± t h < 
^2,h • • ■ < A m> /i and the corresponding eigenfunction approximations u^h, U2,h, ■ • • , u mt h by 
the nonconforming finite element method with the lower bound property Xj^ < Xj (J = 
1,2,-.. ,m). 

Before introducing the postprocessing method, we need to define a finer conforming finite 
element space V h HC which can be constructed by using higher order finite element or refining 
current mesh Th such that this space has the higher convergence order (|19j) 

(5.1) inf llu-T^Ha < Ch 2 \ VnGF 1+7 (fi). 

Now let us introduce a postprocessing method to get the upper bound approximations 
with the help of obtained approximations (Ai^, iti^), • • ■ , (X m ^, u m ,h)- 

Algorithm 5.1. Higher order postprocessing method: 

(1) For j = 1,2, ••■ ,m 

Find Uj.h £ V^ c such that 

(5.2) a(uj, h ,v h ) = Xj ih b( u j,h,Vh), W h eV^ c . 

(2) Construct the finite dimensional space Vh = spanjui^, • ■ ■ , u m ,h} and solve the fol- 
lowing eigenvalue problem in the space V^: 

Find (Aft, u^) ElZ x Vh such that 

(5.3) a(u hj v h ) = X h b(u h ,v h ), \/v e V h . 

Finally, we obtain the new eigenpair approximations (Xi^u^h), - - - , (\ m ,h) ^m,h)- 
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Similarly to the last section, in order to analyze the error estimate for the eigenfunction 
approximations • • ■ , u m> h, we define the projection operator corresponding to the space 
V^ c as follows 

(5.4) a(P* c u,v h ) = a(u,v h ), Vv h e V h HC . 
For this projection operator, we have the following error estimate 

(5.5) \\P^ c u-u\\ a < C inf ||u - v h \\ a < Ch 2 \ 

Theorem 5.2. For the eigenpair approximations (Ai^, u^h), • • ■ , (Kn,h,u m) h) obtained by 
Algorithm \5.1[ we have following estimates 

(5.6) O^Aj-ft-Aj- < Ch*\ 

(5.7) \\u hh -u\\ a < Ch 2 \ 

Proof. First, we prove the error estimate of the eigenfunction approximations (15.71) . From 
(12.21) . (15 .41) and the coercivity of a(-, ■), we have 

\\uj,h ~ Ph C Uj \\l = a(uj, h - Ph C uj, u jth - Ph C Uj) = h[\ h ,,a jh - XjUj, Uj, h - P^Uj) 

— \\^j,h u j,h ~ ^j U j\\o,h\\ u j,h ~ Ph° u j\\a 

< (]>>j,h\\\ U j ~ U j,h\\b + W,h ~ ^j\\\ u j\\b)\\Uj,h ~ Ph° U j\\a 

< Ch^\\u hh -P^ c u 3 \\ a . 
Combined with (1 5 . 5 j) . the following estimate holds 

(5.8) \\%,h-Uj\\a < \\u jth -P? c Uj \\ a +\\PF c u-u\\ a <Cti**. 

Based on the theory in [3] for the error estimate of the eigenvalue problem by finite element 
method, we have the following inequality 

(5.9) - Uj\\ a < C inf \\uj - v h \\ a < C\\uj - u jjh \\ a < Ch 21 . 

v h eV h 

This is the desired result 05 .7p . Now, we come to prove the error estimate property (15. 6p . 
With the help of (14.141) . we have 

(5.10) \\h-\\ < C\\u jth - Uj\\l < Ch 4y . 
Thanks to the minimum-maximum principle ( \2Ah . we have 

(5.11) A,- = min maxi?(t;), 

J V 3 cV,d™Vj=j veVj 

and the discrete version in the space Vh 

(5.12) \j h = min max R(v). 

V Jth cV h ,d™V 3 , h =jv£V J:h 
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Based on (15. lip and f !5. 1 2j) . we can easily obtain Aj < Xj ; h (j = 1, 2, ■ ■ • , m) and complete 
the proof. □ 

6. Numerical Results 

In this section, two numerical examples are presented to validate our theoretical results 
stated in the above sections. 

6.1. Eigenvalue problem on the unit square domain. In this subsection, we solve the 
eigenvalue problem (12. ip on the unit domain Q = (0, 1) x (0, 1). The aim here is to find the 
approximations of the first 6 eigenvalues Ai, A 2 , • • • , A 6 . 

First, ECR element is applied to solve the eigenvalue problem and then the linear finite 
element to do the postprocessing on the series of meshes which are produced by Delaunay 
scheme. The quadratic element is applied to implement Algorithm 15.11 Table [1] shows the 
eigenvalue approximations of the first 6 eigenvalues. And the approximations by postpro- 
cessing method with linear element is presented in Table [2j Table [3] shows the numerical 
results of the postprocessing Algorithm 15.11 with quadratic element. From Table [TJ we can 
find the numerical approximations of ECR element are lower bounds of the exact eigen- 
values. Tables [2] and [3] show the upper bounds of the numerical approximations by the 
postprocessing method using linear and quadratic elements. 

Table 1 . ECR element for eigenvalue problem on unit square 



h 


Ai,/i 




A3, h 




^5,h 




0.2 


19.282902 


46.704835 


46.785588 


72.047988 


88.938276 


89.192096 


0.1 


19.609902 


48.619503 


48.628462 


77.120773 


95.649531 


95.991202 


0.05 


19.711840 


49.180465 


49.181471 


78.517381 


98.019716 


98.037715 


0.025 


19.732395 


49.306198 


49.306523 


78.849396 


98.532500 


98.532943 


0.0125 


19.737526 


49.337705 


49.337756 


78.930339 


98.655173 


98.655315 


Trend 


/ 


/ 


/ 


/ 


/ 


/ 



Table 2. Linear element for postprocessing method on unit square 



h 


Ai,/,, 


^2,h 


A3, h 




^5,h 




0.2 


20.458912 


53.558232 


53.776851 


90.030803 


117.15806 


116.15773 


0.1 


19.940425 


50.568044 


50.561690 


82.204934 


103.84878 


103.61419 


0.05 


19.784700 


49.634441 


49.633984 


79.701583 


99.844154 


99.842216 


0.025 


19.750665 


49.419822 


49.419553 


79.140201 


98.978970 


98.978102 


0.0125 


19.742054 


49.365577 


49.365546 


79.002394 


98.765631 


98.765878 


Trend 


\ 


\ 


\ 


\ 


\ 


\ 
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Table 3. Quadratic element for postprocessing method on unit square 



h 






A3, h 




^5,h 




0.2 


19.744698 


49.427336 


49.434146 


79.315293 


99.263888 


99.326679 


0.1 


19.739654 


49.354899 


49.354991 


78.982884 


98.747751 


98.748129 


0.05 


19.739231 


49.348360 


49.348371 


78.958190 


98.698705 


98.698827 


0.025 


19.739210 


49.348043 


49.348043 


78.956922 


98.696216 


98.696219 


0.0125 


19.739209 


49.348023 


49.348023 


78.956841 


98.696054 


98.696054 


Trend 


\ 


\ 


\ 


\ 


\ 


\ 



Then, EQ™^ element is applied to solve the eigenvalue problem and then the bilinear finite 
element to do the postprocessing on the series of uniform rectangle meshes. Biquadratic 
element is employed to implement Algorithm 15.11 Table H] shows the eigenvalue approxi- 
mations of the first 6 eigenvalues and the approximations by postprocessing method with 
bilinear element is presented in Table [5j Table [6] shows the numerical results of the postpro- 
cessing Algorithm 15.11 with biquadratic element. From Table IU we can find the numerical 
approximations of EQ T ° l element are lower bounds of the exact eigenvalues. Tables [5] and 
[6] show the upper bounds of the numerical approximations by the postprocessing method 
using bilinear and biquadratic elements. 

Table 4. EQ r ° l element for eigenvalue problem on unit square 



h 




^2,h 


^3,h 




^5,h 




1/8 x 1/8 


19.530807 


47.509983 


47.523306 


75.754849 


89.869065 


89.943840 


1/16 x 1/16 


19.686551 


48.878512 


48.879400 


78.123216 


96.394505 


96.399928 


1/32 x 1/32 


19.726009 


49.229994 


49.230051 


78.746202 


98.114256 


98.114609 


1/64 x 1 /64 


19.735907 


49.318474 


49.318478 


78.904035 


98.550189 


98.550211 


1/128 x 1/128 


19.738383 


49.340632 


49.340633 


78.943626 


98.659555 


98.659556 


Trend 


/ 


/ 


/ 


/ 


/ 


/ 



Table 5. Bilinear element for postprocessing method on unit square 



h 




^2,h 


^3,h 




^5,h 




1/8 x 1/8 


20.506335 


52.821611 


54.413322 


90.985945 


114.01626 


114.01577 


1/16 x 1/16 


19.929848 


50.211622 


50.588305 


82.000375 


102.46652 


102.46652 


1/32 x 1/32 


19.786796 


49.563576 


49.656364 


79.717934 


99.633341 


99.633341 


1/64 x 1 /64 


19.751101 


49.401888 


49.424998 


79.147095 


98.930015 


98.930015 


1/128 x 1/128 


19.742182 


49.361487 


49.367259 


79.004399 


98.754514 


98.754514 


Trend 


\ 


\ 


\ 


\ 


\ 


\ 
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Table 6. Biquadratic element for postprocessing method on unit square 



h 




^2,h 


^3,h 




^5,h 


^6,h 


1/8 x 1/8 


19.743647 


49.388394 


49.422024 


79.218966 


99.081741 


99.083312 


1/16 x 1/16 


19.739492 


49.350652 


49.352826 


78.974583 


98.721403 


98.721409 


1/32 x 1/32 


19.739227 


49.348188 


49.348325 


78.957968 


98.697654 


98.697654 


1/64 x 1/64 


19.739210 


49.348032 


49.348041 


78.956906 


98.696145 


98.696145 


1/128 x 1/128 


19.739209 


49.348023 


49.348023 


78.956840 


98.696050 


98.696050 


Trend 


\ 


\ 


\ 


\ 


\ 


\ 



Figure [T] shows the errors of the eigenvalue approximations by ECR and EQ^ elements, 
postprocessing methods with lowest order (linear and bilinear) and higher order (quadratic 
and biquadratic) elements on the unit square. Since we know the exact eigenvalues on the 
unit square, we can give the exact errors of Xj^, A^ and A^ (j = 1,2,3,4,5,6). Since 
the eigenfunctions are smooth, the postprocessing with higher order element can improve 
the convergence order. From Figure [U we can find the eigenvalue approximations have the 
reasonable convergence order. 



10 r- 
10' - 
10° - 
10 1 



10 - 



Errors of eigenvalue on unit square domain 



Errors of eigenvalue on unit square domain 




10 10 1 

Number of elements 



10" 10 10 

Number of elements 



Figure 1. The errors for the eigenvalue approximations on unit square 
by ECR (left) and EQ\ ot (right), where Errl = Ej 6 =i(Aj - A iift ), Err2 = 
ELi(Aj> - Aj) and Err3 = ELlC^ ~ Aj) 



6.2. Eigenvalue problem on the L shape domain. In this subsection, we solve the 
eigenvalue problem (12.11) on the L shape domain Q = (—1, 1) x (—1, 1) \ (—1,0) x (—1,0). 
The aim here is also to find the approximations of the first 6 eigenvalues Ai, A2, • • • , A6- 

First, ECR element is applied to solve the eigenvalue problem and then the linear finite 
element to do the postprocessing on the series of meshes which are produced by Delaunay 
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scheme. The quadratic element is applied to implement Algorithm 15.11 Table [7] shows the 
eigenvalue approximations of the first 6 eigenvalues and the approximations by postprocess- 
ing method with linear element is presented in Table |HJ Table M shows the numerical results 
of the postprocessing Algorithm 15. II with quadratic element. From Table [71 we can find the 
numerical approximations of ECR element are lower bounds of the exact eigenvalues. Ta- 
bles [8] and [9] show the upper bounds of the numerical approximations by the postprocessing 
method using linear and quadratic elements. 

Table 7. ECR element for eigenvalue problem on L shape domain 



h 




^2,h 


^3,h 




^5,h 


^6,h 


0.2 


8.9126839 


14.379736 


18.200634 


26.360803 


27.160361 


34.333785 


0.1 


9.3655553 


14.954062 


19.323660 


28.615430 


30.449112 


39.363462 


0.05 


9.5439611 


15.133230 


19.625778 


29.266600 


31.451422 


40.860929 


0.025 


9.6066551 


15.181407 


19.711040 


29.459256 


31.775515 


41.305926 


0.0125 


9.6277308 


15.193312 


19.732422 


29.506604 


31.869752 


41.426233 


Trend 


/* 


/ 




/ 




/ 



Table 8. Linear element for postprocessing method on L shape domain 



h 


Ai,/, 




A3, h 




A5,/i 


A6,/i 


0.2 


10.578106 


16.642707 


22.460586 


36.034300 


39.441562 


52.785280 


0.1 


9.9724152 


15.610430 


20.469954 


31.132236 


34.109328 


44.593560 


0.05 


9.7438532 


15.309456 


19.934653 


29.958076 


32.522745 


42.389559 


0.025 


9.6728756 


15.224811 


19.786815 


29.627051 


32.086619 


41.719036 


0.0125 


9.6521963 


15.203927 


19.750613 


29.546869 


31.965675 


41.541732 


Trend 


\ 


\ 


\ 


\ 


\ 


\ 



Table 9. Quadratic element for postprocessing method on L shape domain 



h 


Ai,ft 


^2,h 


A3, h 




A5,/i 


A6,/i 


0.2 


9.7067976 


15.233894 


19.813091 


29.809049 


32.402601 


42.472832 


0.1 


9.6691745 


15.201608 


19.744920 


29.541711 


32.011248 


41.586610 


0.05 


9.6496097 


15.197646 


19.739617 


29.522960 


31.939132 


41.497609 


0.025 


9.6432648 


15.197293 


19.739232 


29.521572 


31.921467 


41.481392 


0.0125 


9.6414840 


15.197258 


19.739210 


29.521488 


31.916952 


41.477779 


Trend 


\ 


\ 


\ 


\ 


\ 


\ 



Then, EQ™ 1 element is applied to solve the eigenvalue problem and then the bilinear finite 
element to do the postprocessing on the series of uniform rectangle meshes. Biquadratic 
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element is employed to implement Algorithm 15.11 Table HD1 shows the eigenvalue approxi- 
mations of the first 6 eigenvalues and the approximations by postprocessing method with 
bilinear element is presented in Table [TTJ Table [12] shows the numerical results of the 
postprocessing Algorithm 15.11 with biquadratic element. From Table [10| we can find the nu- 
merical approximations of EQ™* element are lower bounds of the exact eigenvalues. Tables 
[TT1 and [12] show the upper bounds of the numerical approximations by the postprocessing 
method using bilinear and biquadratic elements. 

Table 10. EQ r ° l element for eigenvalue problem on L shape domain 



h 


Al,fc 




A3, h 








1/4 x 1/4 


9.2784846 


15.049120 


19.477978 


28.869068 


30.381898 


39.579205 


1/8 x 1/8 


9.5063501 


15.154836 


19.675337 


29.348136 


31.415853 


40.855247 


1/16 x 1/16 


9.5896364 


15.185968 


19.723326 


29.477224 


31.747401 


41.282853 


1/32 x 1/32 


9.6205870 


15.194336 


19.735243 


29.510335 


31.855244 


41.413996 


1/64 x 1 /64 


9.6323169 


15.196509 


19.738218 


29.518686 


31.891900 


41.454533 


Trend 


/ 


/ 


/ 


/ 


/ 


/ 



Table 11. Bilinear element for postprocessing method on L shape domain 



h 




^2,h 


A3, h 


Azl,h 






1/4 x 1/4 


10.164089 


15.980053 


20.773284 


32.476652 


35.807091 


48.127411 


1/8 x 1/8 


9.7907347 


15.392383 


19.994161 


30.245451 


32.945969 


43.311982 


1/16 x 1/16 


9.6867567 


15.246072 


19.802707 


29.701314 


32.192812 


41.959311 


1/32 x 1/32 


9.6552886 


15.209476 


19.755068 


29.566371 


31.992010 


41.603765 


1/64 x 1 /64 


9.6451377 


15.200312 


19.743173 


29.532700 


31.936225 


41.509772 


Trend 


\ 


\ 


\ 


\ 


\ 


\ 



Table 12. Biquadratic element for postprocessing method on L shape domain 



h 


Ai,/i 


^2,h 


A3, h 


A/l,h 




A6,/i 


1/4 x 1/4 


9.6733499 


15.208409 


19.749420 


29.581107 


32.072149 


41.784684 


1/8 x 1/8 


9.6525125 


15.198129 


19.739857 


29.525433 


31.949976 


41.516052 


1/16 x 1/16 


9.6447652 


15.197334 


19.739249 


29.521742 


31.925488 


41.485165 


1/32 x 1/32 


9.6417221 


15.197261 


19.739211 


29.521499 


31.917575 


41.478309 


1/64 x 1 /64 


9.6405166 


15.197253 


19.739209 


29.521483 


31.914580 


41.475981 


Trend 


\ 


\ 


\ 


\ 


\ 


\ 



Figure [2] shows the errors of the eigenvalue approximations by ECR and EQ^ ele- 
ments, postprocessing methods with lowest order (linear and bilinear) and higher order 
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(quadratic and biquadratic) elements on the L shape domain. Since we don't know the 
exact eigenvalues on the L shape domain, we can only give the errors of A^ — Xj t h and 
Aj /j — Xj h (j = 1,2,3,4,5,6). Since the eigenfunctions here are singular, the convergence 
order by postprocessing with higher order element can not be improved which is shown in 
Figure EJ 



Errors of eigenvalue on L shape domain Errors of eigenvalue on L shape domain 




10 1 10 2 10 3 10 4 10 5 10 1 10 2 10 3 10 4 10 5 



Number of elements Number of elements 

FIGURE 2. The errors for the eigenvalue approximations on L shape domain 
by ECR (left) and EQ\ ot (right), where Errl = £?=i&,* - \ j>h ) and Err2 = 

Yfj=l(^j,h ~ ^j,h) 



7. Concluding remarks 

In this paper, we analyzed the lower bound approximation of eigenvalue problem by 
nonconforming elements (ECR and EQ™*) and also two postprocessing methods to obtain 
the upper bound of the eigenvalues. Especially, based on the lower bound approximations, 
a new postprocessing method which can produce not only higher order convergence but also 
upper bound approximation of the eigenvalues is proposed. This improves the efficiency 
of solving eigenvalue problems and obtain the accurate a posteriori error estimates by the 
lower and upper bounds of eigenvalues. 

We should point out that all the methods and results here can be easily extended to the 
three dimension case. We listed some related space and results. The corresponding ECR 
element in TZ 3 is defined as 

Vf^° : = |f E L 2 {Q) : v\k € span{l, x, y, z, x 2 + y 2 + z 2 }, J v l^ds = j ' v\ K2 ds, 
(7.1) when K x n K 2 = £, and J v\ K ds = 0, if K n dQ = 
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where K, Ki,K 2 G %■ 

The corresponding EQ 1 ^ element in 1Z 3 is defined as 

Vr =={* . : .|« . span {1 , jf, Ms = J^s, 

(7.2) if ^n^ 2 = «, and Jv\ K ds = 0, ifA' nffl = ^}, 

where K,Ki,K 2 G 7&. 

These two nonconforming elements can be used in the three dimensional case to get 
the lower bounds of eigenvalues and the corresponding postprocessing methods can also be 
constructed to obtain upper bounds of the eigenvalues. 

References 

M. Armentano and R. Dufan, Asymptotic lower bounds for eigenvalues by nonconforming finite 
element methods, ETNA 17(2004), 93-101. 

I. Babuska and J. Osborn, Finite element Galerkin approximation of the eigenvalues and eigenvectors 
of selfadjoint problem, Math. Comp., 52(1989), 275-297. 

I. Babuska and J. Osborn, Eigenvalue Problems, in Handbook of Numerical Analysis, V. II: Finite 
Element Methods (Part I), Edited by P. G. Ciarlct and J. L. Lions, 1991, Elsevier. 
P. Ciarlct, The Finite Element Method for Elliptic Problems, North-Holland, 1978. 
F. Chatelin, Spectral Approzimations of Linear Operators, Academic Press, New York, 1983. 
K. Feng, A difference scheme based on variational principle, Appl. Math and Comp. Math, 2(1965), 
238-262. 

J. Hu, Y. Huang and H. Shen, The lower approximation of eigenvalue by lumped mass finite element 
methods, J. Comput. Math., 22(2004), 545-556. 

J. Hu, Y. Huang and Q. Lin, The analysis of the lower approximation of eigenvalues by nonconforming 
elements, to appear, 2010. 

Y. Li, Lower approximation of eigenvalue by the nonconforming finite element method, Math. Numer. 
Sin., 30(2)(2008), 195-200. 

Q. Lin, H. Huang and Z. Li, New expansions of numerical eigenvalues for —Aw = Xpu by noncon- 
forming elements, Math. Comput., 77(2008), 2061-2084. 

Q. Lin and J. Lin, Finite Element Methods: Accuracy and Improvements, Science Press, Beijing, 
2006. 

Q. Lin, L. Tobiska and A. Zhou, On the superconvergence of nonconforming low order finite elements 
applied to the Poisson equation, IMA. J. Numer. Anal., 25(2005), 160-181. 

Q. Lin, H. Xie, F. Luo, Y. Li and Y. Yang, Stokes eigenvalue approximation from below with non- 
conforming mixed finite element methods, Math, in Practice and theory, 19(2010), 157-168. 
Q. Lin, H. Xie, and J. Xu, Lower bounds of the discretization for piecewise polynomials, 



http://arxiv.org/abs/1106.4395l 2011. 

H. Liu and L. Liu, Expansion and extrapolation of the eigenvalue on Q™* element, Journal of Hebei 
University, 23(2005), 11-15. 

H. Liu and N. Yan, Four finite element solutions and comparisions of problem for the Poisson equation 
eigenvalue, J. Numer. Method & Comput Appl, 2(2005), 81-91. 



LOWER AND UPPER BOUNDS OF EIGENVALUE PROBLEM 



19 



[17] R. Rannacher, Nonconforming finite element methods for eigenvalue problems in linear plate theory, 

Numer Math 33, 23-42, 1979 
[18] G. Stang, G. Fix, An Analysis of the Finite Element Method, Englewood Cliffs, NJ: Prentice-Hall, 

1973. 

[19] J. Xu and A. Zhou, A two-grid discretization scheme for eigenvalue problems, Math. Comput., 
70(233)(2001), 17-25. 

[20] Y. Yang, Finite Element Methods Analysis to Eigenvalue Problem, Guizhou People Press, China, 
2004. 

[21] Y. Yang, Z. Zhang and F. Lin, Eigenvalue approximation from below using nonforming finite elements, 
Sci. China Math., 53(1)(2010), 137-150. 

* LSEC, ICMSEC, Academy of Mathematics and Systems Science, Chinese Academy of 
Sciences, Beijing 100190, China 

E-mail address: fshluo@lsec.cc.ac.cn 

t LSEC, ICMSEC, Academy of Mathematics and Systems Science, Chinese Academy of 
Sciences, Beijing 100190, China 
E-mail address: linq@lsec.cc.ac.cn 

* LSEC, ICMSEC, Academy of Mathematics and Systems Science, Chinese Academy of 
Sciences, Beijing 100190, China 

E-mail address: hhxie@lsec.cc.ac.cn 



